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Harvard  Ocean  Prediction  System  (HOPS) 
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Professor  Allan  R.  Robinson 
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Gordon  McKay  Professor  of  Geophysical  Fluid  Dynamics 
Division  of  Engineering  and  Applied  Sciences 
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Pierce  Hall,  29  Oxford  Street 
Cambridge,  MA  02138 

This  is  the  final  technical  report  for  the  project  “Development  of  a  Regional  Coastal 
and  Open  Ocean  Forecast  System:  Harvard  Ocean  Prediction  System  (HOPS)”.  This 
project  was  originally  entitled  “Development  of  a  PE  GULFCAST  Scheme  and  the  Impact 
of  Real  Time  lES  Data  on  GULFCASTS”  and  was  retitled  “Development  of  a  Regional 
Coastal  and  Open  Ocean  Forecast  System”.  The  initial  focus  of  the  project  was  the 
development  of  a  forecast  scheme  based  on  primitive  equation  dynamics  for  the  Gulf  Stream 
region.  During  the  reporting  period,  the  focus  evolved  away  from  a  single  identified  region 
and  a  specific  dynamical  model  towards  a  portable  forecasting  system,  with  the  potential 
to  include  differing  dynamical  models,  for  use  in  any  coastal,  shelf  or  deep  ocean  region. 
The  lead  scientists  at  Harvard  have  been  Dr.  Hernan  Arango  and  Dr.  Patrick  J.  Haley, 
Jr.  with  significant  additional  support  from  Dr.  Carlos  J.  Lozano  and  Wayne  G.  Leshe. 

The  research  tasks  for  this  project  have  been  to:  develop,  calibrate  and  validate 
the  primitive  equation  based  dynamical  model  component  of  the  system;  develop  and 
implement  the  software  for  a  modular,  integrated  forecast  system;  identify  and  acquire 
models  from  the  scientific  community  to  be  utilized  as  modules  within  the  forecast  system; 
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distribute  and  support  the  use  of  HOPS  outside  of  Harvard;  and,  demonstrate  and  validate 
HOPS  via  real  time  exercises  on  sea  and  land. 

General  software  improvements  to  the  HOPS  system  include;  the  addition  of  a  prim¬ 
itive  equation  dynamical  model  to  the  dynamical  model  option;  improvement  in  the  abil¬ 
ity  to  track  software  bug  fixes  and  improvements  (SCCS);  utihzation  of  a  standard  in¬ 
put/output  data  format  (netCDF)  to  increase  product  portability;  and,  implementation 
of  a  procedure  (GNUmake)  to  ease  software  installation. 

Specific  improvements  to  the  primitive  equation  dynamical  model  are  in  the  areas  of: 
ability  to  work  with  high/steep  bottom  topography  (2(7,  Vp);  momentum  flux  computation 
for  terrain-following  coordinates;  vertical  sub-grid  scale  parameterization  (Pacanowski  and 
Philander,  surface  layer).  The  system  is  capable  of  being  utilized  in  coastal  regions  now 
as  coastal  masking  has  been  implemented.  Spherical  coordinates  can  be  used,  as  well  as 
’’rotated  spherical”  coordinates.  Nested  physical  domains  of  varying  horizontal  resolution 
can  share  information  through  n-level,  2-way  nesting.  New  capabilities  include:  addi¬ 
tional  radiation  boundary  conditions  (Orlanski,  Modified  Orlanski);  01  data  assimilation; 
lagrangian  drifter  simulators;  surface  flux  forcing;  and,  coupled  biological  models. 

The  Harvard  Ocean  Prediction  System  (HOPS)  is  now  a  flexible,  portable  and  generic 
system  for  nowcasting,  forecasting  and  simulations.  Recent  reviews  by  Robinson  [1996], 
Robinson  et  al.  [1996b]  and  Lozano  et  al.  [1996]  provide  a  comprehensive  presentation  of 
the  system  and  overview  of  its  applications,  including  coupled  physical-acoustical  [Robin¬ 
son  and  Lee,  1994]  and  physical-biological  [McGillicuddy  et  al.,  1995  a,b]  studies.  A 
rigorous  quantitative  verification  has  been  achieved  for  the  Iceland-Faeroe  Islands  frontal 
system  [Miller  et  al,  1995  a,b;  Robinson  et  al,  1996a].  An  additional  significant  cali¬ 
bration  and  verification  study  of  the  models,  methodology  and  procedures  is  presented  in 
Gangopadhyay  et  al,  (1997);  Robinson  and  Gangopadhyay,  (1997);  Gangopadhyay  and 
Robinson,  (1997).  Important  aspects  of  HOPS  have  been  developed  as  6.1  research  and 
subsequently  transitioned  to  this  6.2  project. 

The  overall  system  schematic  is  shown  on  Fig.  1.  The  heart  of  the  system  for  most 
coastal  applications  is  a  primitive  equation  physical  dynamical  model  which  has  been  spe- 
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cially  structured  for  accurate  and  efficient  calculations  over  steep  topography  [Lozano, 
et  al,  1994;  Haley,  1996].  Vertical  coordinate  options  include  sigma,  hybrid  and  multiple 
sigma  coordinate  transformations  which  are  calibrated  for  specific  applications  via  sensitiv¬ 
ity  analyses  to  both  vertical  and  horizontal  resolutions  [Sloan,  1996,  Chapter  II].  Horizontal 
coordinate  options  include  multiple  two-way  nests  [Sloan,  ibid.].  A  variety  of  physical,  bi¬ 
ological  and  acoustical,  in  situ  and  remotely  sensed,  data  types  have  been  assimilated  in 
a  variety  of  applications.  The  data  analysis  and  management  modules  of  HOPS  represent 
a  major  resource  of  the  Harvard  system.  HOPS  methodology  involves  the  construction  of 
a  best  possible  initial  synoptic  realization  as  a  starting  point  for  the  assimilation  of  new 
synoptic  data.  There  is  an  emphasis  on  the  treatment  of  data  prior  to  assimilation  in  or¬ 
der  to  maximize  the  impact  of  new  data  in  the  light  of  prior  data.  Structured  data  models 
[Lozano  et  al,  1996]  are  utilized  for  this  purpose,  e.g.,  feature  models  or  typical  synoptic 
structures  [Gangopadhyay  et  al..,  1997;  Robinson  and  Gangopadhyay,  1997;  Gangopadhyay 
and  Robinson,  1997]  and  empirical  orthogonal  functions  (EOFs)  in  one  to  three  dimen¬ 
sions.  Recent  developments  include:  i)  the  use  of  temperature  and  salinity  based,  rather 
than  velocity  based,  feature  models  which  are  more  suitable  over  steep  topography;  ii)  the 
addition  of  a  dynamically  balanced  vertical  velocity  to  the  feature  models  for  initialization 
and  assimilations  which  will  couple  physics  and  biology;  and,  iii)  the  optical  dynamical 
component  (Fig.  1).  For  the  construction  of  best  possible  synoptic  realizations  multiple 
data  streams,  including  structured  data  models  and  seasonally  adjusted  historical  synoptic 
realizations,  are  melded.  The  modularity  of  HOPS  facilitates  the  selection  of  a  subset  of 
modules  to  form  an  efficient  configuration  for  specific  applications  and  also  facilitates  the 
addition  of  new  or  substitute  modules. 

A  robust  (suboptimal)  optimal  interpolation  (01)  data  assimilation  scheme  with 
weights  set  by  simple  engineering-type  assumptions  has  been  used  in  HOPS  for  several 
years  [Lozano  et  al.,  1996].  01  was  adopted  in  order  to  focus  research  resources  on  the  as¬ 
similation  of  real  ocean  data  sets  rapidly  into  the  ocean  dynamical  models  for  a  first  round 
of  impact  studies.  A  second  quasioptimal  assimilation  scheme  option.  Error  Subspace  Sta¬ 
tistical  Estimation  (ESSE),  has  been  developed  as  part  of  a  Ph.D.  thesis,  supported  under 
6.1,  and  recently  added  to  the  system  under  the  continuing  6.2  project.  A  rational  ap- 
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proach  was  used  to  identify  an  efficient  statistical  estimation  scheme  feasible  for  use  in  real 
time  with  real  oceanic  data  sets.  The  ESSE  goal  is  to  determine  the  nonlineax  evolution 
of  the  oceanic  state  by  minimizing  the  most  energetic  errors  under  the  constraints  of  the 
dynamical  and  measurement  models  and  both  of  their  uncertainties  [Lermusiaux,  1997]. 
Error  propagation  is  estimated  via  an  ensemble  forecast  using  the  full  nonlinear  model. 
The  evolving  error  subspace  is  characterized  by  singular  error  vectors  and  values,  i.e.,  time 
evolving  three  dimensional  error  EOFs.  Melding  weights  for  assimilation  are  determined 
using  a  minimum  error  variance  criterion.  Importantly,  melding  occurs  in  the  error  sub¬ 
space  and  is  thus  much  less  costly  than  a  classical  analysis  with  the  full  error  covariances. 
The  error  subspace  is  updated  at  the  melding  step  by  combining  the  forecast  principal 
errors,  i.e.,  errors  arising  from  the  dynamical  model  and  the  loss  of  predictability,  with  the 
error  covariances  of  the  measurements. 

Real  time  and  at  sea  forecasts  have  been  carried  out  for  more  than  a  decade  in  twenty 
sites  in  the  Atlantic  and  Pacific  Oceans  and  the  Mediterranean  Sea  (Fig.  2).  Those 
experiments  which  took  place  during  the  time  period  covered  by  this  report  were  conducted 
in:  the  Iceland- Faeroe  Frontal  Region  (1992  and  1993),  the  eastern  Mediterranean  (POEM- 
BC  1995  LIW  Experiment),  the  Strait  of  Sicily  (1994,  1995  and  1996  [NATO  exercise  Rapid 
Response  97]),  the  Haro  Straits  (1996)  and  the  Tyrrhenian  Sea  Skerki  Bank  region  (1996). 
Informal  reports  describing  these  experiments  have  been  prepared  and  are  available. 

HOPS  has  now  been  developed  to  the  level  of  a  true  operational  forecast  system, 
allowing  for  rapid  set-up,  implementation,  and  execution  in  a  new  observational  or  opera¬ 
tional  location  with  Harvard  researchers  functioning  as  a  operational  research  group.  The 
HOPS  system  has  also  been  utilized  by  researchers  at  the  following  locations  (alphabet¬ 
ically):  CNR,  Bologna,  Italy;  IMB,  Ancona,  Italy;  IMS-METU,  Erdemli,  Turkey;  lOLR, 
Haifa,  Israel;  JPL/NASA,  Pasadena;  NCMR,  Athens,  Greece;  NOVA  University;  NRL- 
Stennis;  Naval  Postgraduate  School;  ORI,  Tokyo,  Japan;  Rutgers  University;  SACLANT 
Undersea  Research  Centre,  La  Spezia,  Italy;  Scripps  Institution;  Univ.  of  Massachusetts, 
Dartmouth;  Univ.  of  Rhode  Island;  and.  Woods  Hole  Oceanographic  Institution. 
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